Mon. Not. R. Astron. Soc. 000, [j-^ (1997) Printed 1 February 2008 (MN WT^ style file vl.4) 



Chromatic and spectroscopic signatures of 
microlensing events as a tool for the gravitational 
imaging of stars 

David Valls-Gabaud^'^ 

^ URA CNRS 1280, Observatoire de Strasbourg, 11, rue de I'Universite, 67000 Strasbourg, France 
^ Royal Greenwich Observatory, Madingley Road, Cambridge CB3 OEZ, UK. 



ABSTRACT 

The detection of microlensing events from stars in the Large MageUanic Cloud 
and in the Galactic bulge raise important constraints on the distribution of 
dark matter and on galactic structure, although some events may be due to a 
new type of intrinsic variability. When lenses are relatively close to the sources, 
we predict that chromatic and spectroscopic effects are likely to appear for 
a significant fraction of the microlensing events. These effects are due to the 
differential amplification of the limb and the centre of the stellar disc, and 
present a systematic dependence with wavelength and time that provide an 
unambiguous signature of a microlensing event (as opposed to a new type of 
intrinsic stellar variability). In addition, their measure would provide a direct 
constraint on stellar atmospheres, allowing a 3-dimensional reconstruction or 
imaging of its structure, a unique tool to test the current models of stellar 
atmospheres. 
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1 INTRODUCTION 

Following the pioneering idea by Paczinsky (1986), the detections of several microlensing 
events from stars in the Large Magellanic Cloud (Alcock et al. 1993, 1996; Aubourg et al. 
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1993) and of more than 100 from stars in the Galactic bulge (Udalski et al. 1993, 1994; 
Alcock et al. 1995, 1997) have raised important constraints on the structure of the different 
components of our Galaxy. The associated databases containing millions of light curves are 
also a goldmine for studies in stellar variability (see Ferlet et al. 1997, and references therein). 

Yet doubts may reasonably be cast on the detections of at least some microlensing 
events. The two EROS candidates are variable stars : EROS 1 is a Be star (Beaulieu et al. 
1995), while EROS 2 is an eclipsing binary (Ansari et al. 1995). Although spectra of other 
candidates taken after the event (Delia Valle 1994) or during the event itself (Benetti et al. 
1995; Lennon et al. 1996) show no apparent signs of anomalies, the possibility remains that 
these events are just detections of a new type of variable star. The detection of "bumpers" 
(Alcock et al. 1996), a hitherto unknown type of variable star, and the possibility that 
dwarf novae and cataclismic variables mimick some events (Delia Valle & Livio 1996) are 
particularly worrying, even though the distribution of the events in the HR diagram indicates 
that not all the events can be explained by such phenomena. In addition to continuous 
and follow-up monitoring, it seems important to find unambiguous proofs that the events 
detected can indeed be associated with microlensing events. 

In this context it is interesting to note that the large rate observed towards the bulge 
may be accounted for if about 50 to 90% of the events are due to lenses within the bulge 
itself (Kiraga & Paczyhski 1994), particularly if the bar is oriented close to the line of sight. 
The same applies if many of the lenses are within the LMC (Wu 1994; Sahu 1994) or in the 
halo and bulge of M31. 

If this is the case, a significant number of microlensing events should present character- 
istic chromatic (Valls-Gabaud 1994; Witt 1995; Bogdanov & Cherepaschuk 1995; Gould & 
Welch 1996) and spectroscopic signatures (Valls-Gabaud 1994; Loeb & Sasselov 1995), and 
we present in the Letter detailed predictions for both effects. These effects arise from the 
differential amplification between the limb and the inner regions of the stellar disc, where the 
emerging intensities are different due to the temperature gradient in the atmosphere. The 
differential amplification also produces a non-zero degree of polarisation during the events 
(Simmons et al. 1995ab), and a spectroscopic line shift (Maoz & Gould 1994). 
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2 MICROLENSING OF STARS AS EXTENDED SOURCES 

The amplification produced by a lens of mass Ml at a distance Dql on a source assumed to be 
point-like at a distance Dqs (Refsdal 1964; Liebes 1964) is Ap{y) = {y^ + 2){yy/y'^ + 4)"^ 
where y is the dimensionless distance between the projected position of the lens and the 
centre of the source, y = ^Ja^ + (t — toY/r'^ and ao = Oo/Oe , r = DqlOevJ^ are re- 
spectively the minimum dimensionless impact parameter at time to, and time scale of 
the event, with v± the transverse velocity of the lens. 6'^; is the angular Einstein radius 
6e = {^GMlDi^s/c^DolDos)^^'^ where Dls is the distance separating the lens and the 
source. Since all scales involved are much larger than the wavelength of light, the geometri- 
cal optics approximation is valid, and the events should not show diffraction, interference or 
chromatic effects. Yet, given the colour-magnitude diagrams of the bulge and the LMC, and 
the sensibility limits of the present surveys, one expects that between one and two thirds 
of the events are associated with giant stars, so that the point source approximation might 
be questionable. If the angular radius 9^ of a star of radius R^^ seen at a distance Dqs is 
comparable with the angular Einstein radius 6e of the lens, then effects arising from the 
finite size of the stars will become apparent. Thus it is only when either the lens mass Ml 
or the relative distance Di^s/Dqs is small (or both) that these angles are similar. Although 
the conditions for this to be true might seem contrived, there will always be a number of 
microlensing events where the lens will transit through the stellar disc (Gould 1994; Witt 
1995). The fraction of all lensing events which have a maximum impact parameter smaller 
than the Einstein radius is simply J-'=<6^><1/6-e>, where average must be taken over 
the distribution of distances of both lenses and sources, and hence depends on the galactic 

— 1/2 

model adopted. Since the fraction increases like M^ , if the lenses are sub-stellar objects 
the fraction may reach unity for a large range of Dls distances. However for lenses of 0.1 
Mq and typical fractional distances Di^s/Dqs from 0.1 (bulge) to about 5 10~^ (M31), the 
fraction of transit events goes from 31% to 16% for typical clump giants of 10 Rq. Since they 
constitute at least a third of the sources, the minimum fraction of transit events is around 
10% for the bulge and about 5% in M31, for 0.1 lenses. For these events, the sizes of 
the stars has to be taken into account. 
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Figure 1. Chromatic signatures of the microlensing of a giant red star of TcH = 4000 K, logg=2.0 and solar metaUicity, for 
the indicated geometrical parameters [70 = 7(fo)]. The magnifications for the extended Johnson photometric broad bands are 
calculated assuming a square root non-linear limb darkening profile. The lower panels give the relative magnification relative 
to the one in the Johnson V band, Q = Ax/ Ay which refiects the relative change in colour index A(m;\ — mv) ~ 1.086(1 — Q). 
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3 CHROMATIC AND SPECTROSCOPIC SIGNATURES 

The inclusion of the finite size of the sources was first analysed by Bontz (1979) and more 
recently by Nemiroff & Wickramasinghe (1994), Gould (1994) and Witt & Mao (1994) for 
the idealised case of a uniform brightness source, but little attention has been paid so far to 
the amplification produced by actual stars. Yet it is well known that even in the Eddington 
approximation the difference in intensity between the centre and the limb of a star reaches 
almost a magnitude. In the general case then, the magnification A^, of an extended source 
with specific intensity X(tz7) is given by the integration over the extension f2 of the source 
of the amplification produced on an infinitesimal area on the source at an angular distance 
9 from the point lens Ae{0) = dw X{w) Ap{6) / dzu X{w). Assuming that the 
star is spherically symmetric {i.e., neglecting rotational fiattening, spots, etc.) and using 
polar coordinates in the source plane and the notation = 7^ + + 2 7 x cos the actual 
magnification of the star becomes 



where x = r/R^, a* = O^jO-^, 7 = 9i^os/9^ is the angular separation between the lens and 
the centre of the star in units of the stellar angular radius, Bx{x) = Xa(x)/Xa(0) is the 
normalised brightness profile, and Rx = 2 J dx x B\ is the dimensionless effective radius at 
that wavelength. In the case where Bx{x) is constant, the maximum amplification reduces 
to A"^^^^^^ = (l + 4/a2)i/2 (^j^efs^ial 1964, Bontz 1979), while for 7 > 1, ^ Ap, the 
point-source limit. The brightness profile Bx is given by the solution of the radiative transfer 
equation : 



where Xa(0) is the emerging intensity at the centre of the star and Sx the source function. It 
is the explicit dependence with wavelength of the brightness profile Bx and the differential 
amplification that produce both chromatic and spectroscopic effects (Valls-Gabaud 1994), 
that reflect the temperature gradient T{tx) in the stellar atmosphere. 

Given a stellar atmosphere model, one derives the corresponding limb proflles for both the 
continuum and the lines. We have used the large grid of LTE models from Kurucz (1994) and 
fltted, for the both continuum and a selected set of lines, simple analytical expressions. The 
reason for this is twofold : flrst the computational speed is substantially increased, and second 
the information content of the limb darkening is limited and prevents the use of very detailed 
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models (Bohm 1961). Since the pioneering studies of Milne (1922), the limb profiles for the 
continuum are approximated by a linear law i3^™(x) = 1 — xa (1 — /u), where /i = (1 — x^)^/^. 
In the modelling of stellar atmospheres, however, non-linear laws have yielded significantly 
better fits and in particular the logarithmic law BY~"^{x) = 1 — xa (1 — /i) — yxfiln^ and 
the square root law (Diaz-Cordoves et al. 1995) B^^^{x) = 1 — xx{l — fi) — y\ (j- — y/J^j 
seem to provide good fits to most of the computed profiles (Van Hamme 1993), although we 
stress that actual precise data is only available for the Sun. In order to predict the different 
light curves in the photometric passbands used in microlensing monitoring projects, the 
limb darkening coefficients are integrated over these bands, and our results agree with those 
published by Claret & Gimenez (1992) and Van Hamme (1993), for both the Johnson and 
Stromgren systems. 

The important point to note is that since the limb darkening coefficients obviously depend 
on wavelength due to the variation in opacity, different limb profiles will appear depending 
on the photometric band. This strong dependence with wavelength implies that the same 
microlensing event will present quite different light curves in different photometric bands, 
as shown in Fig. |^. 

We can understand the behaviour of the light curves as a function of wavelength by noting 
that the effective radius R\ increases from 0.64 in the U band to 0.90 in the K band for the 
typical red giant star in Fig. |l[ Thus, the longer the wavelength, the more uniform the star 
is, while in the near UV stars are more sharply peaked, giving a larger amplification when 
the impact parameter is small. At large impact parameters the strong decrease in brightness 
at shorter wavelengths takes over and the magnification becomes smaller, relative to larger 
wavelengths. Using the actual limb darkening profiles, the linear Milne law or the non-linear 
profiles do not make any significant difference {AA/A < 0.01) in the final magnifications, 
and hence the use of the linear approximation is fully justified. 

One should stress that all these predictions are based on theoretical models which have 
never been tested in such detail. There may well be further complications, such as for 
example, at larger wavelengths, in the far IR or sub-mm range, where the limb can actually 
be brighter (at least in the Sun). In extended envelopes, the increased opacity in the UV can 
produce a larger Rx in this domain, and decrease the contrast Q, while in the X-ray domain, 
the coronal emission is very patchy and a smooth, spherically symmetric distribution is 
meaningless. There may be also complications related to the presence of chromospheres, 
which change the temperature gradient and hence the limb profiles. The spectral lines will 
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Figure 2. Spectroscopic signatures of microlensing, for the same giant star considered in Fig. and two typical lines : H/3 
and CO 2.3 /xm. On the right panels the effect of increasing the stellar diameter (from a* = 1 to a* = 2) is illustrated for the 
otherwise same transit conditions. 

be more sensitive to this effet (see below). To first order, the chromatic effect does not 
depend on the colour of the stellar source, since the integrated intensity of the star does 
not appear in Eq. |l|. However, there will be some dependence since limb darkening profiles 
do depend on the spectral type of the star, as well as on the luminosity class. In general, 
the limb darkening increases when the gravity or temperature decreases and also when the 
metallicity increases, so the giants of the Galactic and M31 bulges are ideal probes. 

The best way to measure the chromatic effect is to observe the event in photometric 
passbands as widely separated in wavelength as possible, ideally U and K which are the 
most discriminatory for giant stars. It is unfortunate that the current surveys use bands as 
close as V , R and /, where the effect is predicted to be very small. As Fig. |] shows, the 
relative magnifications Q (with respect to the V band) are smaller than 10%, that is, the 
change in colour index A(mA — my) ~ 1.086(1 — Q) is less than 0.1 magnitudes, hence it 
is also essential to do differential, rather than absolute, photometry to measure the effect. 



© 1997 RAS, MNRAS 000, [[jl^ 



8 D. Valls-Gahaud 

Present technology may achieve 0.001 mag (Young et al. 1991), although a precision of 0.01 
will be enough. Finally, the use of narrow-band Stromgren photometry is also recommended, 
since the Stromgren bands probe different opacity regimes where the source function changes 
strongly. The wavelength and time dependence of the light curves are unlike those produced 
by any intrinsic stellar variability, and hence constitute a unique signature of microlensing 
events. 

The spectroscopic signature arises, once again, because line profiles change significantly 
across the disc, although in this case the variation will depend on the mechanism of formation 
of the line. For example strong resonant lines, like Ly-a or Ca II H+K, Mg II h+k present 
a decreasing relative depression with respect to the continuum, and hence appear brighter 
at the limb, while pure absorption lines disappear at the limb. Lines are far more sensitive 
to the actual gradient of temperature in the atmosphere, since they probe a larger range 
of depths, while others are more sensitive to pressure or turbulence. However, since the 
observations of the variability of the line profiles seems difficult with the present technology, 
a more appropriate observable is the equivalent width (EW) of the line. For very strong lines, 
narrow-band photometry in the appropriate range could also be used (Loeb & Sasselov 1995), 
but the availability of real-time spectra makes this less practical. To illustrate the predicted 
behaviour of the EWs during a microlensing event[j]. Fig. |^ presents the predictions for the 
Balmer H/9 and CO 2.3 /im lines, the CO forming at lower temperature, much higher in the 
atmosphere, and is limb-brightened. It is immediately apparent that the amplitude of the 
spectroscopic effect may reach 20%. In the case of H/?, which is limb-darkened, the EW first 
decreases when most of the fiux originates from the limb, then increases when the lens is 
probing the central areas of the star. For the vibrational-rotational bands of CO (from 2.2 to 
4.6 /im), the opposite behaviour is present, and the temperature gradient is probed at higher 
altitude in the atmosphere. Once again, the temporal variation of the lines during the event 
constitute a unique signature of microlensing, unlike any variation produced by intrinsic 
stellar variability. Note also that the time scale is no longer the Einstein ring crossing time, 
but rather the stellar crossing time ~ R^v^{l — Di^s/Dqs) (see Fig. |^). 

t Detailed predictions for any type of star of metallicity between 0.002 and 0.04 in the wavelength range 1200A-5/im are 
available upon request. 
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4 DISCUSSION 

The gravitational imaging of stars is then possible, where not only the brightness distribution 
on the stellar disc is measured (via the solution of the integral equation in Eq. |1|), but also 
the source function can be recovered, solving Eq. 0, hence giving a 3-dimensional picture of 
the stellar atmosphere. This is unprecedented in astrophysics, where only the Sun and very 
few nearby stars are spatially resolved, either directly (e.g., Gilliland & Dupree 1996) or by 
interferometry (e.g., Baldwin et al. 1996). 

Actual techniques for the inversion of Eqs. |l| and |^ cannot fit within the scope of this 
Letter, and will be dealt with elsewhere. We simply note here the important implications 
of stellar imaging. The predictions presented here are based on model atmospheres, in LTE 
or NLTE, and the expected limb darkening coefficients have seldom been compared with 
observations others than the solar ones. The reason is that the only way to measure these 
coefficients (besides the solar case) accurately is using eclipsing binaries, which must have 
circular orbits, well-detached components, and produce total eclipses. Even then, the derived 
coefficients correlate with the assumed sizes of the stars (e.g., Tabachnik 1969), so the 
measure of limb darkening coefficients in microlensing events is a unique opportunity to test 
stellar atmosphere models. Note however that the inclusion of non-linear terms produces 
insignificant differences in the resulting magnification, yet even the measure of the linear 
coefficients is important, for they provide an estimation of the contribution function of the 
continuum (basically the integrand of Eq. This in turn gives a direct estimation of the 
temperature profile, something that is known in very few stars (e.g. Matthews et al. 1996) 
and that provides unprecedented constraints on the opacity sources. This technique, applied 
to the lines (ideally the profiles, but the EWs are still extremely important) may lead to the 
element abundance imaging of the star. In general, atmospheric diagnostics will be possible, 
by selecting lines sensitive to pressure, or temperature, etc. 

Unlike Doppler or Zeeman-Doppler imaging, the gravitational imaging is unbiased to- 
wards bright stars with large spots, although we note that the presence of spots will consid- 
erably complicate the inversion unless the periods are much longer than the time scale of the 
event. Complications may also arise from blending (Kamionkowsky 1995) and extinction. 
The ongoing real-time follow-up teams, PLANET (Albrow et al. 1997) and GMAN (Abe 
et al. 1997), could adapt their observing strategies to measure these signatures, by increas- 
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ing the time sampling, widening the wavelength range covered and obtaining the highest 
resolution spectra possible. 

In conclusion, the chromatic and spectroscopic signatures of microlensing events allow 
us not only to deduce properties of the lenses, but also, and maybe more importantly, to 
open an entire new window on stellar atmospheres. 
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Note: After this paper was completed, Alcock et al. (1997, ApJ, submitted, astro-ph 9702199) have reported real-time pho- 
tometry and spectroscopy of events associated with a M4 giant in the Galactic bulge. Both signatures may be present, although 
further analysis is required for confirmation. 
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